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Foreword 


This  report  was  prepared  as  part  of  the  Statistical  Methods  for  Drug  Testing  project  (Program 
Element  0305889N,  Work  Unit  030.‘5889N.R2143DR001).  sponsored  by  the  Chief  of  Naval 
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metliodologies  for  the  analysis  of  drug  testing  programs  and  data.  The  work  described  here  was 
performed  during  FY93. 

Tltis  is  the  third  in  a  series  of  reports  on  the  use  of  Markov  chains  lor  the  analysis  of  random 
urinalysis  programs.  The  first  report  is  Markov  Chains  for  Random  Urinalysis  I:  Age-Test  Model 
(NPRDC-TN'93-5).  The  second  report  is  Markov  Chains  for  Random  Urinalysis  11:  Age-Test 
Model  With  Absorbing  State  (NPRDC-TN-93-6).  Related  work  also  includes  Probability  of 
Detection  of  Drug  Users  by  Random  Urinalysis  in  the  U.S.  Navy  (NPRDC-TN-93-2). 

We  thank  Cecil  Hombeck  from  the  San  Diego  Navy  Dmg  Screening  Laboratory  for  the  use  of 
his  data. 


MURRAY  W.  ROWE 
Director,  Manpower  Systems  Department 


v 


Summary 


Background 

Since  198 1 ,  the  U.S.  Navy  has  maintained  a  zero  tolerance  drug  policy.  The  cornerstone  of  this 
policy  is  an  aggressive  urinalysis  testing  program.  The  testing  program  is  intended  to  deter  and 
detect  drug  abuse,  as  well  as  provide  data  on  the  prevalence  of  drug  abuse. 

Previous  work  (Thompson,  Boyle,  &  Hentschel.  1993;  Boyle,  Hentschel,  &  Thompson,  1993) 
developed  Markov  models  for  analyzing  random  urinalysis  strategies  stratified  by  “time  last 
tested"  (i.e.,  age-test  strategies).  This  report  describes  a  method  to  quantify  the  effects  of  daily 
variations  in  random  urinalysis  testing  rates  on  detection  of  di'ug  users.  The  report,  an  extension  of 
Thompson  and  Boyle  (1993),  presents  a  daily  Markov  model  that  includes  drug  excretion  rate 
kinetics. 

Objective 

The  objective  of  this  work  is  to  quantify  the  effects  of  daily  variations  in  random  urinalysis 
testing  rates  on  detection  of  drug  users. 

Approach 

First,  daily  drug  testing  is  modeled  as  a  Markov  chain.  This  chain  includes  parameters  for  daily 
drug  testing  probabilities  and  for  conditional  probabilities  of  testing  positive.  Second,  a  model  for 
the  conditional  probabilities  of  testing  positive  for  cocaine  use  is  presented.  Drug  dosage  is  treated 
as  a  random  variable  and  its  distribution  is  developed  using  Navy  data.  Third,  times  to  detection 
for  various  diug  testing  and  diug  usage  patterns  are  compared  using  the  above  models. 

Results 

Analyses  of  four  examples  of  weekly  testing  patterns  are  presented.  The  first  two  examples 
represent  the  obser\'ed  random  sampling  behavior  of  two  Navy  activities.  The  third  example  shows 
equal  testing  rates  on  weekdays  and  no  testing  on  weekends.  The  fourth  example  is  equal  testing 
rates  every  day  of  the  week.  All  of  the  examples  have  been  standardized  to  a  15%  monthly  testing 
rate. 

In  each  example,  once  per  week  cocaine  use  is  assumed.  The  average  time  to  detection  for  the 
first  activity  varies  from  12.7  to  409.7  months  depending  on  the  day  of  the  week  cocaine  is  used. 
For  the  second  activity  ,  average  time  to  detection  varies  from  12.2  to  106.5  months.  If  an  activity 
does  not  test  on  weekends,  a  Friday  night  user  averages  271.1  months  to  detection.  Average  time 
to  detection  under  constant  rate  daily  testing  is  27.3  months. 

Conclusions 

A  Markov  model  for  random  drug  urinalysis  testing  that  allows  for  daily  variations  in  testing 
probabilities  was  developed.  The  formulation  allows  for  any  fixed  length  cycle  (e.g.,  week,  month). 
Drug  kinetics  and  drug  user  gaming  can  be  incorporated  into  the  Markov  model  via  conditional 


probabilities.  The  Markov  chain  provides  estimates  of  the  distribution  of  time  to  detection  and 
mean  time  to  detection. 

The  cocaine  kinetics  model  presented  here  provides  tlie  probability  of  testing  positive  as  a 
function  of  time  since  ingestion.  Initial  dose  is  treated  as  a  random  variable  and  tlie  distribution  of 
initial  dose  is  estimated  from  Navy  data. 

The  analyses  have  shown  that  time  to  detection  varies  dramatically  with  varying  (observed) 
daily  testing  rates.  Unequal  daily  testing  rates  provide  opponunities  for  gaming  drug  users  to 
extend  the  mean  time  to  detection.  Gaming  is  not  possible  with  equal  probabilities  of  testing  across 
days. 
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Introduction 


Since  1981 ,  the  U.S.  Navy  has  maintained  a  zero  t->i mce  drag  policy.  ITie  cornerstone  of  this 
policy  is  an  aggressive  urinalysis  testing  program.  ring  program  is  intended  to  deter  and 

detect  drug  abuse,  as  well  as  provide  data  on  the  presale'  e  of  drug  abuse.  All  officer  and  enlisted 
personnel  are  subject  to  random  urinalysis  testing  on  a  continuing  basis.  Current  policy  (Chief  of 
Naval  Operations  [CNO],  1990)  directs  Navy  commands  to  test  10  to  20%  of  their  members  each 
month.  The  Navy’s  random  urinalysis  program  has  been  successful  in  reducing  the  proportion  of 
service  members  testing  positive  for  drugs  (7%  in  1983  compared  to  1%  in  1991).  However,  since 
any  drug  abuse  impacts  readiness,  health,  and  safety,  continuous  evaluation  and  improvement  of 
the  Navy’s  program  is  required. 

Previous  reports  (Thompson.  Boyle,  &  Hentschel.  1993;  Boyle.  Hentschel,  &  Thompson  1993) 
developed  Markov  models  for  analyzing  random  urinalysis  strategies  stratified  by  "time  last 
tested’’  (i.e.,  age-test  strategics).  In  related  work,  Evanovich  (1985)  uses  Markov  chains  to  model 
the  probability  of  detection  of  drug  users.  Stoloff  (1985)  investigates  the  relationship  between 
deterrence  and  random  urinalysis  testing. 

This  report  describes  a  method  to  quantify  the  effects  of  daily  variations  in  random  urinalysis 
testing  rates  on  detection  of  drug  users.  The  repon,  an  extension  of  Thompson  and  Boyle  (1992), 
presents  a  daily  Markov  model  that  includes  drug  excretion  rate  kinetics.  The  Markov  models 
allow  for  a  fixed  length  daily  cycle  (e.g.,  v/eekly,  monthly)  and  include  an  absorbing  state  for 
detection.  Dnig-specific  urinary  excretion  rates  are  also  modeled.  Initial  drug  dose  is  treated  as  a 
random  variaole.  Optimal  gaming  strategies,  using  data  from  the  current  Navy  urinalysis  program, 
are  then  developed.  These  strategies  treat  initial  drug  dose  as  a  random  variable.  Finally,  the  report 
presents  the  survival  distributions  of  time  to  detection  under  optimal  gaming  and  under  constant 
testing  rates. 

A  relevant  issue,  not  previously  addressed,  is  drug  detection  time  in  urine.  TTie  detection  times 
vary  by  drug.  For  example,  marijuana  is  detectable  for  approximately  2  to  4  days,  cocaine  for 
approximately  1  to  2  days,  and  methamphetamine  for  approximately  2  to  3  days.  Detection  times 
also  vary  by  dose,  analytical  test  method  used,  individuals  physical  condition,  fluid  intake,  and 
method  and  frequency  of  ingestion.  Discussions  of  urinary  excretion  kinetics  of  cocaine, 
marijuana,  methamphetamines,  and  amphetamines  appear  in  Ambre,  Ruo,  Nelson,  and  Belknap 
(1988);  Beckett  and  Rowland  (1965);  Hamilton,  et  al.  (1977);  Johansson,  Gillespie,  and  Halldin 
(1990);  Johansson  and  Halldin  (1989);  and  Cook,  et  al,  (1992). 

Existence  of  drug  user  gaming  can  be  gleaned  from  the  Department  of  Defense  surveys  of 
substance  abuse  and  health  behavior  among  militaiy'  personnel.  Five  surveys  have  been  conducted 
since  1980.  TTie  most  recent  was  the  1992  Worldwide  Survey  (Bray,  et  al.,  1992).  These  surveys 
provide  comprehensive  and  detailed  estimates  of  the  prevalence  of  drug  use  among  military 
personnel.  Results  from  the  1992  Worldwide  Survey  show  that  55.9%  of  drug  users  believe  "drug 
users  cunail  use  when  tliey  think  they  will  be  selected  for  urinalysis.”  About  6.6%  of  surveyed 
Navy  personnel  admitted  to  some  drug  use  in  the  past  12  months.  Tiis  indicates  a  need  for  random 
urinal}  sis  programs  that  are  not  predictable  or  programs  that  drugs  users  do  not  believe  are 
predictable.  As  a  result,  drug  user  guming  sU'ategies  are  investigated  here. 
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The  Daily  Model 


Using  llu*  theory  and  notation  developed  in  previous  repons  (Thompson,  et  al.,  1 993;  Boyle,  et 
al.,1993),  a  class  of  Markov  chains  with  transition  matrices.  P.  of  the  following  form  can  be 
defined: 


0  1  -ai/>|  0  0 

0  0  1  '  a2/>2  0 


0  a,/), 

0  a2P2 


P« 


0 

1  - 

0 


0 

0 

0 


0  0 
0  0 
0  0 


0  ajfj 

0  1 


(1) 


Given  that  an  individual  is  in  state  /,  p,  is  the  probability  of  the  individual  being  tested  and  a,  is  the 
conditional  probability  of  a  positive  result.  Note  that  an  individual  is  in  state  /  (2  <  j  <  cO  if  an 
individual  was  in  state  i  - 1  one  time  period  earlier  and  eitlier  was  not  tested  or  tested  negative.  An 
individual  is  in  state  1  if  an  individual  was  in  state  done  time  period  earlier  and  either  was  not  tested 
or  tested  negative.  State  if  +  1  is  tlie  single  absorbing  state.  An  individual  is  in  slate  d  +  1  (detected) 
if  tested  positive.  For  the  applications  in  this  report,  a  fixed  length  daily  cycle  of  d  -  7  is  assumed. 
Here  states  1  to  7  represent  Sunday  through  Saturday,  respectively. 


The  p/s  and  a,’s  are  allowed  to  be  unrestricted  probabilities.  That  is,  0  <  p,  <  1  and  0  <  U/  <  1 , 
V/.  However,  a^p,  4^  0  for  at  least  one  t.  This  constraint  guarantees  that  states  1  through  d  are 
transient  and.  starting  in  any  transient  state,  an  individual  must  be  caught  in  finite  time  with 
.  ^nainty. 


With  the  above  assumption  and  a  theorem  in  Taylor  and  Karlin  (1984).  the  matrix  P  from 
Equation  1  can  be  svritten  as: 


P- 


Q  R 
O  I 


with 


0  1  -  a^p^  0  0  •  •  •  0 

0  0  1  -  U2P2  0  -  -  •  0 

0  0  0  1  -  «3P3  •  •  •  0 
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1  ■  ^d\Pd\ 
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and  I  -  Q  nonsingular.  Recalling  tliat  the  matrix 


W-(I-Q)' 


is  the  fundamental  matrix,  some  standard  results  needed  in  tlic  next  section  can  be  stated.  First,  tlie 
following  notation  is  employed: 


W,^  is  the  //th  element  of  W. 

T  is  time  to  detection  or  absorption  by  state  J  +  1 , 
Py  is  the  /^h  element  of  P",  and 
is  the  state  of  the  process  at  time  n. 

Using  these  definitions  and  theory  from  Taylor  and  Karlin  (1984): 


7^  <  I  .^0  ”  * 

£[7'|X,  =  ij-tw„-T„ 


1 


and 

{VAR  jr  I  Xo-/j  }-(2W-l)T-tjp 


for  I  -  1 . d.  XsQ  is  the  elementwise  square  of  x.  If  the  row  vector  Hq  is  any  initial  distribution 

for  the  process  concencrated  on  the  transient  states,  the  above  implies 


P 


T<n 


d 


-I  n„(i)P',., 

/=  1 


and  the  unconditional  distribution  function  of  7' can  be  computed.  Its  expectation  and  variance  are 


L, 


-Hox 


I 


and 
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Finally,  for  avi  individual,  Uic  number  of  tests  Y  over  d  time  peiiods  has  mean  and  variance 


E\Y 


I-  I 


VAR 


Y 


I 


1-  1 


r/1  -P,)- 


Drug  Kinetics 


Gibaldi  and  Perrier  (1982)  characterize  drug  elimination  from  the  body  after  rapid  intravenous 
injection  as  following  first-order  kinetics.  Specifically,  the  rate  of  loss  of  the  drug  from  the  body  is 
given  by 


-  -kXit) 
at 


(2) 


where  X  (r)  is  the  amount  of  drug  in  tlie  body  at  time  /  al'tcr  injection.  The  constant  A  is  called  the 
apparent  first  order  elimination  rate.  Equation  2  is  a  first  order  linear  differential  equation  and  is 
easily  solved  to  yield 


Xij)  -  xoe 


wlieic  JTo  anioum  of  dnig  ingested  (i  e..  the  initial  dose). 

Additionally,  some  of  the  drug  is  excreted  unchanged  in  the  urine  (i.e.,  renal  excretion)  and 
again  by  first-order  kinetics  me  rate  of  appearance  of  intact  drug  In  the  urine  is  proixjrtional  to  the 
amount  of  drug  in  (he  body.  1  hus, 

(3) 


where  A'„(;)  is  the  cumulative  amount  of  drug  excreted  in  the  urine  to  time  t  and  A,  is  the  appai'cnt 
first-order  rate  constant  for  renal  excretion.  Integrating  Equation  3  yields 

Note  that  limj^  A'„(r)  <-  k^x^A  and  A/k  represents  the  fraction  of  the  initial  dose  eliininateo 
unchanged  in  the  urine. 

Drug  elimination  from  the  body  tlirough  metabolism  also  follows  first-order  kinetics  with 

=kjX(t)-k„Mir) 
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where  A/(r)  is  ilie  amount  of  metabolite  in  tlic  body  at  time  ?.  Tne  constants  k,  and  are  die 
respective  anpareni  first-order  rate  constants  for  metabolite  formation  and  elimination.  Equation  5 
is  linear  and  by  elementjny  methods  can  be  solved  to  yield 


The  differential  equation  describing  the  appearance  of  metabolite  in  the  urine  is 


where  Mu  (/;  is  the  amount  of  metabolite  in  the  urine  and  k„u  is  the  apparent  first  order  rate  constant 
for  the  excretion  of  metabolite  in  tlie  mine.  Using  Equations  6  and  7  yields  the  time  rate  of  change 
of  metabolite  in  the  urine  as 


dMu 

dt  “  k-kS 


e  ■  *»' .  e-  -  f  f 


Ambre  (1985)  subjected  previously  published  urinary  excretion  data  on  cocaine  and  its 
metabolites  to  kinetic  analysts  using  the  first-order  models  described  above.  He  obtained  the 
following  parameter  estimates  for  cocaine  (COC)  and  its  metabolite  benzoylecgonine  (BZ): 


k  =k 


0.464 

0.03 

0.2134 

0.0923. 


These  values  can  be  used  in  Equations  3  and  8  yielding  time  paths  of  change  in  milligrams  per  hour. 
Ambre  also  assumed  a  urination  rate  of  1  milliliter  per  minute.  Thus,  concentrations  in  nanograms 
per  milliliter  are  obtained  through  multiplication  of  Equations  3  and  8  by  the  factor  lOOO^  'OO.  The 
resulting  expressions  are 


^COc(ri  •“ 


lOOO^  (.03)  (.464) 


-  232 XQe ' 


1000'  (.0923)  (.2134) 


•  092,V  ^  -■464/ 


60  (.464  -  .0923) 


-  883.1863^0  c 


,  -  .OT2  V  _  _  -  .464/ 
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where  and  are  respectively  urine  concentrations  for  COC  and  BZ  in  nanograms  per 

milliliter  (ng  ml)  at  r  hours  after  ingestion  of  a  dose  jr^in  milligrams  (mg)  of  cocaine.  To  illustrate. 
Figure  1  is  a  somiloganihmic  plot  (converted  to  logs  base  10)  of  concentration  curves  for  dose 
levels  of  50  and  100  milligrams  of  cocaine. 


0  12  24  36  48  60  72 


Time  (hours) 

Figure  1.  Concentrations  of  cocaine  and  benaoyiecgonine 
in  urine  versus  titne  from  ingestion. 

Concentrations  of  cocaine  and  benzoylecgonine  for  1.076  urinalysis  tests  were  obtained  from 
the  San  Diego  Navy  Drug  Screening  Labotatory.  This  represents  all  the  specimens  that  tested 
positive  (or  BZ  and  had  detectable  amounts  of  COC  over  a  2-year  period.  Figure  2  is  a  scatter  plot 
of  BZ  concentration  versus  COC  concentration  on  these  paired  values. 


10*  10^  10*  itr*  10*  10* 

Cocaine  Concentration  (ng/ml) 

Figure  2.  Scatter  plot  of  benzoylecgonine  concentration  versus  cocaine 
concentration  in  urine. 


Equations  9  and  10  can  be  solved  for  the  initial  dose  and  the  time  since  ingestion  r  by  taking 
natural  logs  and  simplifying.  The  resulting  two  linear  equations  in  two  unknov.'ns  in  matrix  form 
are 


—  -“I 

— 

1  -.464 

Inxo 

.  in 

Ycoc/232 

1 

1  -.0923 

t 

In 

1(883.1863^000  +  232Tbz)  <  (232  •  883.1863) 

L  J 

Solving  the  above  equations  for  each  of  the  1,076  observations  yields  a  distribution  on  initial  dose 
xq.  Figure  3  is  a  frequency  distribution  of  log,o  jcc  in  milligrams. 
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Figure  3.  Frequency  distribution  of  logio  initial  dose  in  miliigranis. 
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The  probabilities  of  testing  positive  can  now  be  determined  from  the  distribution  of  initial 
dose  .To.  At  any  time  t  after  ingestion,  the  probability  of  testing  positive  is 


a(r) 


P 


1bz(02:  150ng/ml 


#  of  Xq  values  with  ^BZ 


t)>  150 


1076 


The  probabilities  of  testing  positive  for  values  r  -  12, 36, 60 .  156  hours  were  computed.  These 

times  correspond  to  1  to  7  days,  assuming  drug  ingestion  at  10:00  p.m.  and  testing  at  10:00  a.m. 
each  day.  Table  1  lists  the  probabilities. 
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Table  I 


Probabilities  of  Testing  Positive  for  Cocaine  Use 
1  to  7  Days  After  Ingestion 


Day 

1  (houw> 

a(;) 

1 

12 

.9376 

2 

36 

.6508 

3 

60 

.1229 

4 

84 

0 

5 

108 

0 

6 

132 

0 

7 

156 

0 

Applications 


Analyses  of  four  examples  of  weekly  testing  patterns  are  presented  in  tliis  section.  ITie  first  two 
examples  (UIC  A  and  UiC  B)  represent  the  observed  random  sampling  behavior  of  two  Navy 
activities.  The  results  are  averaged  over  a  2'year  period.  The  third  example  shows  eqital  testing 
rates  on  weekdays  and  no  testing  on  v/eekends.  The  fourth  example  is  equad  testing  rates  eveiy  day 
of  the  week.  All  of  the  examples  have  been  standardized  to  a  15%  monthly  testing  rate,  or  a  3.5% 
weekly  testing  rate  (7  /  30  •  0.1 5  -  0.035).  A  monthly  rate  of  i5%  represents  the  midpoint  of  the 
mandated  (CNO,  19^)  Navy  range  of  monthly  rates  (10  to  20%). 


Table  2  contains  the  weekly  random  samplmg  testing  pattern,  the  probability  of  being  tested, 
and  the  mean  time  to  detection  by  day  of  week  for  UIC  A.  The  first  line  of  values  in  the  table 
(Percentage  of  Tests)  is  the  observed  testing  pattern  expressed  as  the  percentage  of  tests  conducted 
by  day  of  week.  The  second  line  (Probability  Tested)  is  the  probability  of  being  tested  standardized 
to  a  0.035  weekly  testing  rate.  The  third  line  (Average  Detection  Time  [months])  is  the  average 
time  to  detection  assuming  cocaine  is  used  once  a  week  and  on  the  same  day  of  the  week  each 
week.  UIC  A  does  almost  all  of  its  testing  on  Saturday,  Sunday,  and  Monday.  If  cocaine  is  used 
every  week  on  Friday  night,  the  average  time  to  detection  is  12.7  months.  If  cocaine  is  used  every 
week  on  Saturday  night  the  average  time  to  detection  is  13.3  months.  Compare  this  to  cocaine  use 
on  Monday  night  where  the  average  time  to  detection  is  409.7  months  or  34  years.  While  the  UIC 
A  strategy  catches  the  weekend  useis,  it  allows  gaming  drug  users  to  significantly  increase  time  to 
detection. 


Table  2 

Daily  Model  for  Random  Sample  Testing  and  Cocaine  Use  (UIC  A) 


Day  of  the  Week 


Sun 

Mon 

l\iC 

Wed 

Thu 

Fri 

Sat 

Percentage  of  Tlrsts 

29.4 

34.5 

0 

2.5 

0 

2.5 

31.1 

Probability  Tested 

.0103 

.0121 

0 

.0009 

0 

.0009 

.0109 

Average  Detection 
Time  (months) 

20.3 

409.7 

251.4 

122.3 

25.4 

12.7 

13.3 

Another  example  of  an  observed  Navy  testing  pattern  is  shown  in  Table  3.  This  activity  does 
half  of  its  testing  on  Fridays.  The  best  strategy  for  a  once  a  week  cocaine  user  at  this  command  is 
to  u‘  e  on  Sunday  with  an  average  detection  time  of  106.5  months  or  8.9  years.  Table  4  shows  the 
effects  of  not  testing  on  week  nds.  Cocaine  could  be  used  once  a  week  on  Friday  with  a  271. 
month  average  time  to  detection.  Table  5  shows  the  average  time  to  detection  when  the  daily  testing 
rate  is  constant  across  day  of  week.  Regardless  of  when  the  drug  user  ingests  cocaine  average  time 
to  detection  is  27.3  months. 


Tables 

Daily  Model  for  Random  Sample  Testing  and  Cocaine  Use  (UiC  B) 


Day  of  the  Week 

Sun 

Mon 

Tiie  Wed 

Thu 

Fri 

Sat 

Percentage  of  Tests 

8.8 

0 

7.1  13.3 

11.5 

50.5 

8.8 

Probability  Tested 

.0031 

0 

.0025  .0047 

.0040 

.0177 

.0031 

Average  Detection 
Time  (months) 

106.5 

39.8 

25.5  14.9 

12.2 

47.7 

73.0 

Table  4 

Daily  Model  for  Random  ^  unple  Testing  and  Cocaine  Use 
(Constant  Testing  on  Weekdays) 

Day  of  tlie  Week 

Sun 

Mon 

TUe  Wed 

Thu 

Fri 

Sat 

Percentage  of  Tests 

0 

20 

20  20 

20 

20 

0 

Probabilir/  Tested 

0 

.007 

.007  .007 

.007 

.007 

0 

Average  Detection 
Time  (months) 

19.4 

19.4 

19.4  20.9 

35.5 

271.1 

43.0 

Table  5 


Daily  Model  for  Random  Sample  Testing  and  Cocaine  Use 
(Constant  Testing  on  Everyday) 


Day  of  tlie  Week 

Sun 

Mon 

Tuc 

Wed 

Thu 

Fri 

Sat 

Percentage  of  Tests 

14.3 

14.3 

14.3 

14.3 

14.3 

14.3 

14.3 

Probability  Tested 

.005 

.005 

.005 

.005 

.005 

.005 

.005 

Average  Detection 

Time  (months) 

27.3 

27.3 

27.3 

27.3 

27.3 

27.3 

27.3 

Conclusions 


A  Markov  model  for  random  drug  urinalysis  testing  that  allows  for  daily  variations  in  testing 
probabilities  was  developed.  For  example,  if  a  command  does  not  test  on  weekends,  this  mod' !  can 
be  used  to  investigate  the  effect.  Although  the  examples  presented  here  are  weekly  models,  tlie 
formulation  allows  for  any  fixed  length  cycle  (e.g.,  month).  Drug  kinetics  can  be  incorporated  into 
the  Markov  model  via  conditional  probabilities.  The  conditional  probabilities  of  testing  positive 
given  that  an  individual  is  tested  are  included  as  model  parameters.  Also,  drug  user  gaming 
strategies  can  be  included  via  these  conditional  probabilities.  The  Markov  chain  provides  estimates 
of  the  distribution  of  time  to  detection  and  mean  time  to  detection. 

Tire  cocaine  kinetics  model  presented  here  provides  the  probability  of  testing  positive  as  a 
function  of  time  since  ingestion.  Initial  dose  is  treated  as  a  random  variable  and  the  distribution  of 
initial  dose  is  estimated  from  Navy  data.  Models  for  other  drugs  could  be  developed  from  the 
literature  cited  earlier. 

The  analyses  have  shown  that  time  to  detection  varies  dramatically  with  varying  (observed) 
daily  testing  rates.  Unequal  daily  testing  rates  provide  opportunities  for  gaming  drug  users  to 
extend  the  mean  time  to  detection.  Gaming  is  not  possible  with  equal  probabilities  of  testing  across 
days. 


10 


References 


Ambre.  J.  (1985).  The  urinary  excretion  of  cocaine  and  metabolites  in  humans:  A  kinetic  analysis  of 
published  data.  Journal  of  Analytic  Toxicology,  9, 241-345. 

Ambre,  J.,  Ruo,  T.  1.,  Nelson,  J.,  &  Belknap,  .S.  ( 1 988).  Lirinar>’  excretion  of  cocaine,  benzoylecgonine, 
and  ecgonine  methyl  ester  in  human.  Journal  of  Analytic  Toxicology,  12, 301-306. 

Beckett,  A.  H.,  &  Rowland,  M.  (1965).  Urinary  excretion  kinetics  of  amphetamine  in  man.  Journal  of 
Pharmacy  and  Pharmacology,  11,  628-639. 

Boyle,  J.  P.,  Hentschel,  D.  J.,  &  Thompson,  T.  J.  (1993).  Markov  chains  for  random  urinalysis  II: 
Age-test  model  with  absorbing  state  (NPRDC-TN-93-6).  San  Diego,  CA;  Navy  Personnel 
Research  and  Development  Center. 

Bray,  R.  M.,  Kroutil,  L.  A.,  Luckey,  J.  W..  Wheeless.  S.  C..  lannacchione,  V.  G.,  Anderson,  D.  W., 
Marsden,  M.  E.,  &  Dunteman,  G.  H.  (195‘2).  1992  worldwide  survey  of  substance  abuse  and 
health  behaviors  among  military  personnel  (R'TT/51 54/06- 16FR).  Research  Triangle,  NC: 
Research  Triangle  Instimte. 

Chief  of  Naval  Operations  (1990).  Alcohol  and  drug  abuse  prevention  and  conttol  (OPNAVINST 
5350.4B):  Author. 

Cook,  C.  E.,  Jeffcoat,  A.  R.,  Sadler,  B.  M.,  Hill,  j.  M.,  Voyksner,  R.  D.,  Pugh,  D.  E..  White.  W.  R„  & 
Perez-Reyes,  M.  (1992).  Pharmacokinetics  of  oral  methamphetamine  and  effects  of  repeated  daily 
dosing  in  humans.  Drug  Metabolism  and  Disposition,  20(6),  856-862. 

Evanovich,  P.  (August  1985).  A  model  for  drug  testing  (CRM  85-33).  Alexandria,  VA:  Center  for 
Naval  Analyses. 

Gibaldi,  M.,  &  Perrier.  D.  (1982).  Pharmacokinetics  (2nd  Ed).  New  York:  Marcel  Dekker,  Inc. 

Harrdlton,  H.  E.,  Wallace,  J.  E.,  Shimek,  E.  L.,  Land,  P.,  Harris.  S.  C.,  &  Christenson,  J.  G.  (1977). 
Cocaine  and  benzoylecgonine  excretion  in  humans.  Journal  of  Forensic  Sciences,  22(4),  697-707. 

Hoel,  P.  G.,  Port,  S.  C.,  &  Stone,  C.  J.  (1972).  Introduction  to  stochastic  processes.  Boston;  Houghton 
Mifflin  Company. 

Johansson,  E.,  Gillespie,  H.  K.,  &.  Halldin,  M.  M.  (1990).  Human  urinary  excretion  profile  after 
smoking  and  oral  administmtion  of  -tetrahydrocannabinol.  Journal  of  Analytical 

Toxicology,  14,  176-180. 

Johansson,  E.,  &  Halldin,  M.  M.  (1989).  Urinary  excretion  half-life  of  A’  -ietrahydrocannabinol-7-oic 
acid  in  heavy  marijuana  users  after  smoking.  Journal  of  Analytical  Toxicology,  13, 218-223. 

Stoloff,  P.  H.  (August  1985).  The  effectiveness  of  urinalysis  as  a  deterrent  to  drug  use  (CNR  111). 
Alexandria,  VA:  Center  for  Naval  Analyses. 


11 


Taylor.  H.  M.,  8r.  Karlin,  S.  (1984).  An  introduction  to  stochastic  modeling.  Orlando,  FL:  Academic 
Press,  Inc. 

lliompson,  T.  J.,  &  Boyle,  J.  P.  (1992).  Probability  of  detection  of  drug  users  by  random  urinalysis 
in  the  U.S.  Navy  (NPRDC-TN-93-2).  San  Diego.  CA:  Navy  Personnel  Reseai'ch  and 
Development  Center. 

Thompson,  T.  J.,  Boyle,  J.  R,  &  Ilentschel.  D.  J.  (1993).  Markov  chains  for  random  urinalysis  I: 
age^test  model  (NPRDC-TN-93-5).  San  Diego,  CA;  Navy  Personnel  Research  and 
Development  Center. 


12 


Distribution  List 


Distribution; 

Chief  of  Naval  Personnel  (PERS-6),  (PERS-63).  (PERS-63E),  (PERS-6C),  (PERS-6E) 
Chief,  Bureau  of  Medicine  and  Surgery,  ^Vashington,  DC  (MED-244) 

ODASD,  Drug  Enforcement  Plans  and  Support 
OCDEP&S,  Science  and  Testing  Programs  Officer 
Defense  Technical  Information  (DTIC)  (4) 

Copy  to: 

Armstrong  Laboratory,  Operations  and  Support  Directorate  (AL/DO),  Brooks  AFB,  TX 
AI7HR-D0KL  Technical  Library,  Brooks  AFB,  TX 
Chief  of  Naval  Operations  (N091 ),  (N8 1 3) 

Chief  of  Naval  Personnel  (PERS-00),  (PERS-OOB).  (PERS-OOW).  (PERS-OIJJ) 

Deputy  Assistant  Secretary  of  the  Navy  (Manpower  and  Reserve  Affairs) 

Director  of  Research,  U.S.  Naval  Academy 
Office  of  the  Assistant  Secretary  of  Defense  (P&R) 

President,  Naval  War  College  (Code  E- 1 1 2 1 ) 

Superintendent,  Naval  Postgraduate  School 


13 


